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Cavity quantum electrodynamics (cavity QED) 
describes the coherent interaction between mat- 
ter and an electromagnetic field confined within 
a resonator structure, and is providing a use- 
ful platform for developing concepts in quan- 
tum information processing^. By using high- 
quality resonators, a strong coupling regime can 
be reached experimentally in which atoms co- 
herently exchange a photon with a single light- 
field mode many times before dissipation sets 
in. This has led to fundamental studies with 
both microwave^'^ and optical resonators^. To 
meet the challenges posed by quantum state 
engineering^ and quantum information process- 
ing, recent experiments have focused on laser 
cooling and trapping of atoms inside an optical 
cavity®''"'^. However, the tremendous degree of 
control over atomic gases achieved with Bose- 
Einstein condensation^ has so far not been used 
for cavity QED. Here we achieve the strong cou- 
pling of a Bose-Einstein condensate to the quan- 
tized field of an ultrahigh-finesse optical cavity 
and present a measurement of its eigenenergy 
spectrum. This is a conceptually new regime of 
cavity QED, in which all atoms occupy a single 
mode of a matter-wave field and couple identi- 
cally to the light field, sharing a single excitation. 
This opens possibilities ranging from quantum 
communication^— liiii^ to a wealth of new phe- 
nomena that can be expected in the many-body 
physics of quantum gases with cavity-mediated 
interactions 

The coherent couphng of a single two-level atom with 
one mode of the quantized light field leads to a splitting 
of the energy eigenstates of the combined system and 
is described by the Jaynes-Cummings modeli^. For the 
experimental realization the strong coupling regime has 
to be reached, where the maximum coupling strength go 
between atom and light field is larger than both the am- 
plitude decay rate of the excited state 7 and that of the 
intracavity field k. In the case of a thermal ensemble of 
atoms coupled to a cavity mode, the individual, position- 
dependent coupling for each atom has to be taken into 
account. 

To capture the physics of a Bose-Einstein condensate 
(BEC) coupled to the quantized field of a cavity, we con- 
sider N atoms occupying a single wave function. Because 
the atoms are in the same motional quantum state, the 
coupling g to the cavity mode is identical for all atoms. 
Moreover, bosonic stimulation into the macroscopically 
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FIG. 1: Experimental situation, a, 36 mm above the cav- 
ity, 3.5 X 10^ ultracold atoms are loaded into the dipole po- 
tential of a vertically oriented one-dimensional optical lat- 
tice. This trumpet-shaped standing wave has its waist inside 
the ultrahigh-finesse cavity and is composed of two counter- 
propagating laser beams. A translation of the lattice trans- 
ports the atoms into the cavity mode. There, they are loaded 
into a crossed-beam dipole trap formed by a focused beam 
oriented along the y axis and one of the transport beams, b. 
Almost pure condensates with 2.2 x 10^ atoms are obtained. 



populated ground state should largely reduce the scat- 
tering of atoms into higher momentum states during the 
coherent evolution. This situation is therefore well de- 
scribed by the Tavis-Cummings modeli^, where N two- 
level atoms are assumed to identically couple to a single 
field mode. A single cavity photon resonantly interact- 
ing with the atoms then leads to a collective coupling of 

A key characteristic of the coupled BEC-cavity system 
is its eigenenergy spectrum, which we map out with a 
single excitation present. An ensemble of thermal atoms 
does not fulfill the requirement of identical coupling, but 
it shows a similar energy spectrum, which can be mod- 
elled by the Tavis-Cummings hamiltonian with an effec- 
tive collective coupling^^. In previous measurement o^^i^^ 



2 



and also in a very recent report^, these eigenenergies 
have been measured for thermal atoms coupled to a cav- 
ity. Aside from the sensitivity of the spectrum to the 
precise spatial distribution of the atoms, the differences 
between a BEC and a thermal cloud, or between a BEC 
and a Mott insulator, should also be accessible through 
the fluctuations of the coupling, that is, in the width of 
the resonances^^. 

First experiments bringing together BEC physics and 
cavities concentrated on correlation measurements using 
single atom counting^^, studied cavity enhanced super- 
radiance of a BEC in a ring cavity^^, observed nonlinear 
and heating effects for ultracold atoms in an ultrahigh- 
finesse cavity2^i2^, and achieved very high control over 
the condensate position within an ultrahigh-finesse cav- 
ity using atom chip technology's. 

To create a BEC inside an ultrahigh-finesse optical cav- 
ity we have modified our previous set-up^. The experi- 
ment uses a magnetic trap 36 mm above the cavity, where 
we prepare 3.5 x 10^ *^Rb atoms in the \F, mp) — |1, —1) 
state with a small condensate fraction present. The 
atoms are then loaded into the dipole potential of a ver- 
tically oriented standing wave, formed by two counter- 
propagating laser beams. By varying the frequency dif- 
ference S between the upwards and the downwards prop- 
agating wave, the standing- wave pattern, and with it the 
confined atoms, move downwards at a velocity v = X 6/2, 
where A is the wavelength of the trapping laser— i^i. Be- 
cause of continuous evaporative cooling during the trans- 
port, the number of atoms arriving in the cavity is re- 
duced to typically 8.4 x 10^ atoms with a small conden- 
sate fraction present. During the 100 ms of transport a 
small magnetic field is applied to provide a quantization 
axis and the sample remains highly spin-polarized in the 
1) state. However, owing to ofF-resonant scattering 
in the transport beams, a small fraction of the atoms 
undergoes transitions into the \F — 2) hyperfine state 
manifold. 

At the position of the cavity mode, the atoms are 
loaded into a crossed-beam dipole trap formed by one of 
the transport beams and an additional, horizontal dipole 
beam with a waist radius of Wx — Wz = 27iim (see 
Fig.[T]). A final stage of evaporative cooling is performed 
by suitably lowering the laser power to final trapping fre- 
quencies {lOx , U, 



= 27r X (290, 43, 277) Hz, ending up 
with an almost pure condensate of 2.2 x 10^ atoms. 

The ultrahigh-finesse cavity has a length of 176 pm and 
consists of symmetric mirrors with a 75 mm radius of 
curvature, resulting in a mode waist radius of 25 pm. A 
slight birefringence splits the resonance frequency of the 
empty cavity for the two orthogonal, principal polariza- 
tion axes by 1.7 MHz. With the relevant cavity QED pa- 
rameters {go, K, j) = 2ttx (10.6, 1.3, 3.0) MHz, the system 
is in the strong coupling regime. The length of the cav- 
ity is actively stabilized using a laser at 830 nm (ref. 26). 
The intracavity intensity of the stabilization light gives 
rise to an additional dipole potential of 2.4 Ercc, with the 
recoil energy defined as i?rcc = h"^ /{2m\^), where m is 
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FIG. 2: Cavity transmission for the cr+ and o~ polar- 
ization component. The probe laser frequency is scanned 
at a speed of 25 MHz ms~^ while the cavity detuning is fixed. 
The original transmission data, recorded with a resolution of 
0.4 ps, is averaged over 2 ms using a sliding average. A single 
peak for each polarization can clearly be distinguished from 
the background of about 60 dark counts per second. 



the mass of the atom. The chemical potential /i — 1.8i?i cc 
of 2.2 X 10^ trapped atoms is comparable to the depth 
of this one-dimensional lattice, so that long-range phase 
coherence is well established in the atomic gas. The 1/e 
lifetime of the atoms in the combined trap was measured 
to be 2.8 s. 

To find the eigenenergies of the coupled BEC-cavity 
system for a single excitation, we perform transmis- 
sion spectroscopy with a weak, linearly polarized probe 
laser of frequency Wp. To this end, the resonance fre- 
quency of the empty cavity is stabilized to a frequency 
LiJc, which in general is detuned by a variable frequency 
Ac ~ LOc ~ LOa, with rcspcct to the frequency of the 
\F ^ I) ^ \F' = 2) transition of the Da line of ^^Rb. 
The transmission of the probe laser through the cavity 
is monitored as a function of its detuning Ap — ujp — uj;^ 
(see Fig. [5]). The two orthogonal circular polarizations 
of the transmitted light are separated and detected with 
single-photon counting modules. The overall detection 
efficiency for an intracavity photon is 5%. To probe the 
system in the weak excitation limit, the probe laser in- 
tensity is adjusted such that the mean intracavity pho- 
ton number is always below the critical photon number 
nn = 7V(25o) = 0-04. A magnetic field of 0.1 G, ori- 
ented parallel (within 10 %) to the cavity axis provides a 
quantization axis. 

From individual recordings of the cavity transmission 
as shown in Fig. [2] we map out the low-excitation spec- 
trum of the coupled system as a function of Ac (see 
Fig. [3]). After resonant excitation we do not detect an in- 
fluence on the BEC in absorption imaging for large atom 
numbers (see Fig.[T]). For small BECs of the order of 
5,000 atoms we observe a loss of 50 % of the atoms after 
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FIG. 3: Energy spectrum of the coupled BEC-cavity 
system. The data points are measured detunings of reso- 
nances for cr^ (red circles) and a~ (black triangles) polarized 
light. Each data point is the average of three measurements 
with an uncertainty of about 25 MHz. The solid lines are the 
result of a theoretical model (see Methods). Bare atomic res- 
onances are shown as dotted lines, whereas the empty cavity 
resonance of the TEMoo mode is plotted as a dashed-dotted 
line. Note the asymmetry in the splitting at Ac = caused 
by the influence of higher-order transverse modes. Neglecting 
this influence, the eigenenergies shown by the dashed lines 
would be expected where the free parameters were adjusted 
to fit the spectrum for Ap < 0. 



resonant probing. The normal mode splitting at Ac = 
amounts to 7 GHz for ct+ polarization, which results in a 
collective cooperativity oi C — Ng^ /{2^k) = 1.6 x 10^. 
The splitting for the cr^ component is smaller, because 
the dipole matrix elements for transitions starting in 
|1, — 1) driven by this polarization are smaller than those 
for (T+. 

A striking feature of the energy spectrum in Fig. [3] is 
a second avoided crossing at probe frequencies resonant 
with the bare atomic transitions \F = 2) \F' = 1, 2, 3). 
It is caused by the presence of atoms in the |F = 2) 
hyperfine ground state. This avoided crossing is lo- 
cated at a cavity detuning where the eigenenergy branch 
of the BEC-cavity system with no atoms in \F = 2) 
would intersect the energy lines of the atomic transitions 
\F — 2) — ^ \F' — 1,2,3). Accordingly, the avoided cross- 
ing is shifted by approximately Ng"^ /Ap = 2tt x 1.8 GHz 
with respect to the intersection of the empty cavity reso- 
nance with the bare atomic transition frequencies. From 
a theoretical analysis (see Methods), we find the size of 
the \F — 2) minority component to be 1.7% of the total 
number of atoms. 

Our near-planar cavity supports higher-order trans- 
verse modes equally spaced by 18.5 GHz, which is of the 
order of the collective coupling g^/N in our system. In 
general, the presence of one additional mode with the 



same coupling but detuned from the TEMqo mode by At 
would shift the resonance frequencies at A^. to first order 
by iV(7^/(2At). In our system, this results in a clearly 
visible change of the energy spectrum with respect to a 
system with a single cavity mode only (see dashed lines in 
Fig. [3]). This can be seen as a variant of the "superstrong 
coupling regime" in which the coupling between atoms 
and the light field is of the order of the free spectral range 
of the cavity. 

We describe the BEC-cavity system in a fully quan- 
tized theoretical model (see Methods), which yields the 
eigenenergies of the coupled system. Good agreement 
between the measured data and the model is found (see 
Fig. [3]) for 154,000 atoms in the |1,— 1) state and for 
2,700 atoms distributed over the Zeeman sublevels of the 
\F = 2) state, with the majority in |2,— 1). The sub- 
stantial influence of the higher-order transverse modes is 
modelled for simplicity by the coupling of the BEC to 
one additional effective cavity mode. 
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FIG. 4: Shift of the lower resonance of the coupled 
BEC-cavity system from the bare atomic resonance. 

The cavity was locked at Ac = 0. and a~ polarization are 
shown as red circles and black triangles, respectively. Each 
data point is the average of three measurements. The atom 
number was determined separately from absorption images 
with an assumed error of ±10%; the vertical error bars are 
too small to be resolved. The dashed lines are fits of the 
square root dependence on the atom number, as predicted 
by the Tavis-Cummings model. The solid lines are fits of a 
more detailed theoretical model (see Methods) resulting in 
maximum coupling rates = 2n x (14.4 ± 0.3) MHz and 
= 27r X (11.3 ± 0.2) MHz. The ratio of the two coupfing 
rates of 1.27 ± 0.03 agrees with the expected ratio of 1.29 of 
the corresponding Clebsch-Gordan coefficients. 

To test the square-root dependence of the normal 
mode splitting on the number of atoms in the BEC, a 
second measurement was conducted. We set the cav- 
ity frequency to Ac — and record the detuning of 
the lower coupled state from the bare atomic resonance 
\F — 1) \F' = 2) as a function of the number of 
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atoms as displayed in Fig.[4l The atom number was var- 
ied between 2,500 and 200,000, determined from sepa- 
rately taken absorption images with an estimated statis- 
tical error of ±10%; possible systematic shifts are esti- 
mated to be within ±7 %. The dependence of | Ap| on the 
number of atoms is well described by a square root, as ex- 
pected from the Tavis-Cummings model (dashed lines). 
However, for a weakly interacting BEC the size of the 
atomic cloud — and thus the spatial overlap with the cav- 
ity mode — depends on the atom number. Our more de- 
tailed model, which includes this effect, as well as the 
influence of higher-order cavity modes, yields maximum 
single-atom couplings of —2ttx (14.4 ±0.3) MHz and 
g„- = 2ttx (11.3±0.2) MHz for the two polarization com- 
ponents (solid lines). The ratio of these two couplings is 
1.27 ± 0.03 and agrees with the ratio of 1.29 that is ob- 
tained from the effective Clebsch-Gordan coefficients for 
the and cr^ transitions starting in state |1, — 1). 

The coupling of a single mode of a matter- wave field to 
a single cavity mode opens a route to new experiments. 
It facilitates the manipulation and study of statistical 
properties of quantum-degenerate matter-wave fields 
by a quantized optical field, or even the generation 
of entanglement between these two field o^^'^^ . The 
detection of single atoms falling through the cavity 
has already been demonstrated with this set-up^. In 
principle, the detection of small impurity components 
embedded in a large BEC presented here can also be 
extended to single atoms. This is an important step 
towards the realization of schemes aiming at the cooling 
of qubits immersed in a large BEG^. 

METHODS SUMMARY 

Optical transport. The transport of the atoms into 
the cavity is accomplished in T = 100 ms with a maxi- 
mum acceleration of a = 22.4 ms^^. The standing wave 
used to transport the atoms has its waist {wx,Wy) — 
(25, 50) pm centred inside the cavity and is locked to a 
i33Qg resonance at a wavelength of 852 nm. The inten- 
sity and frequency of each beam are precisely controlled 
by acousto-optical modulators, which are driven by two 
phase-locked, homebuilt direct digital synthesis genera- 
tors. The frequency difference 6 between the two counter- 
propagating waves follows (5 (t) — [1 — cos (27ri/T)](5max/2, 
with a maximum detuning of ^max = 1,670 kHz. With 
the maximally available power of 76 mW per beam the 
trap depth at the position of the magnetic trap is 1.1 pK. 
During transport, the power in the laser beams is kept 
constant until the intensity at the position of the atoms 
has increased by a factor of ten. Subsequently, this in- 
tensity is kept constant. 

Theoretical model. To gain understanding of the 
presented measurements we have developed a fully quan- 
tized theoretical description of the coupled BEC-cavity 
system. Our model includes all Zeeman sublevels in the 
5^Si/2 and 5^P3/2 state manifolds of ^^Rb and both 
orthogonal polarizations of the TEMqo cavity mode. 
For simplicity, the effect of higher-order transverse 



cavity modes is modelled by the effective coupling to 
one additional mode which is detuned from the TEMqo 
mode by At = 2tt x 18.5 GHz. The free parameters 
of the model are the coupling strength between this 
effective mode and the BEC and the population of 
the several ground states in the condensate. Good 
agreement with the measured energy spectrum is found 
for the ground-state population given in the main text 
and a coupling between BEC and effective cavity mode 
that is r = 1.2 times the coupling to the TEMqo mode. 

METHODS 

To characterize the coupled BEC-cavity system theoreti- 
cally, we start with the second-quantized hamiltonian de- 
scribing the matter-light interaction in the electric-dipole 
and rotating- wave approximation^^. We take all hyper- 
fine states including their Zeeman sublevels in the 5^Si/2 
and 5^P3/2 state manifolds of ®^Rb into account and de- 
scribe the cavity degrees of freedom by two orthogonal 
linear polarizations of the TEMqo mode. We choose the 
quantization axis, experimentally provided by a small 
magnetic field, to be oriented parallel to the cavity axis. 
The near-planar cavity supports higher-order transverse 
modes equally spaced by At = 27r x 18.5 GHz. To incor- 
porate the coupling to all higher-order transverse modes 
we include one additional effective cavity mode in our 
model with its resonance frequency shifted by At with 
respect to that of the TEMqo mode. 

Considering only one spatial atomic mode for the 
ground state manifold and another one for the excited- 
state manifold, the hamiltonian of the uncoupled system 
reads 

1 

i j fc=Qp=^,t 

where the indices i and j label the states 5'^Si/2\F,mp) 
and 5'^F^/2\F' ,mF'), respectively. The operators gj (or 
gi) and (or ej) create (or annihilate) an atom in the 
mode of the corresponding ground and excited states 
with frequencies uig- and cuej ■ The operators o|. ^ (or ak,p) 
create (or annihilate) a photon with energy hiuk and lin- 
ear polarization p in the cavity mode fc, where fc = 0, 1 
labels the TEMqq mode and the additional effective cav- 
ity mode, respectively. 

The coupling between the BEC and the cavity is de- 
scribed by the interaction hamiltonian 

1 

-ffint = X! X! Sif 4 ak,p gi + h.c. 

fe=Op=— «j' 

where g^j;^ denotes the coupling strength for the transi- 
tion i j driven by the cavity mode k with polarization 
p, and h.c. is the hermitian conjugate. For the TEMqq 
mode the coupling strength g^^^ depends on the dipole 
matrix element V^j for the transition i j driven by 
the polarization p, the mode volume Vq, and the overlap 
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Uq between the two spatial atomic modes and the TEMqo 
mode: 
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We numericaUy calculated the ground state of the Gross- 
Pitaevskii equation in the potential formed by the dipole 
trap and the cavity stabilization light for N = 2 x 10^ 
atoms, and found the spatial overlap between BEG and 
TEMqo mode to be Uq = 0.63. Because of the po- 
sition uncertainty of the BEG relative to the cavity 
mode, the overlap might deviate by up to 20% from 
this value. The repulsive interaction between the atoms 
leads to a slight decrease in Uq with increasing atom 
number, which was found numerically to follow Uo{N) — 
VO^il - 0.0017 N°-^^). For the coupling strength to the 
additional effective cavity mode we assume gl^^ 
with r being a free parameter. 

The initial atomic population of the several ground 
states i given by the free parameters Ni is relevant 
for the form of the energy spectrum of the coupled 
BEG-cavity system. To find the energy spectrum in 
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the weak excitation limit we can restrict the analysis 
to states containing a single excitation. These are the 
states jlfe^p; A^i, . . . , iVg; 0), where one photon with po- 
larization p is present in the mode k and all atoms 
are in the ground-state manifold, and also all possible 
states \Ok,p', Ni, . . . , Ni ~ 1, . . . , N^; Ij), where no photon 
is present and one atom was transferred from the ground 
state i to the excited state j. Diagonalization of the 
hamiltonian H = Hq + Hint in the truncated Hilbert 
space spanned by these states yields the eigenspectrum 
of the coupled system as a function of the cavity detuning 
Ac. The relevant transition energies with respect to the 
energy of the initial ground state are plotted in Fig.[3l 

The detuning of the lower resonance branch at Ac = 
is a function of the atom number and given by 



Fitting this dependence to the data in Fig. 3] yields the 
maximum coupling strengths f;^± for the two polarization 
components cr^. 
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